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ABSTRACT 

The details of extended physical processes, such as the gas dynamic 
flow over an airfoil, the reactive flow through a combustor, or the electric 
field in a multi-contact transistor, are understood by solving the dif- 
ferential equations of a mathematical model of the process. The accuracy 
of finite difference methods for the numerical solution of the equations 
is increased if the underlying mesh fits the region boundaries and is 
closely spaced in regions where the solution is rapidly varying. Automatic 
methods for producing a satisfactorily adjusted mesh have been developeed 
for one-dimensional problems. In one simple, effective scheme of this 
kind the unknown function and the distribution of mesh modes are found 
simultaneously, the nodes being placed so that they correspond to points 
uniformly spaced on the solution curve. 

In a two-dimensional generalization, the nodes correspond to points 
equally spaced on the solution surface in two directions that are as nearly 
orthogonal as possible. Examples of such meshes are shown for given surfaces 
in the figures. The meshes fit the circular boundaries and come closer 
together where the given surface is steeper. 
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GENERAL GRID GENERATION METHOD 


Map problem domain onto unit cube in 
computational coordinate space. Uniform, 
rectangular grid in computer coordinates 
gives curvilinear grid in original domain. 



Choose map to reduce truncation error of 
finite difference soiution scheme for 
problem. 


Problem: Find z(x,y) so that 

*xx “*■ 2yy = P^Jt,y) In + y^ < 1. 

z = b{x,y) on x^ + y^ = 1 . 

with p and b given functions. 

Change coordinates: (x,y) — ( ^ , tj ) 

{*) _i_ ~ = io 

9^ j 9r, j 

ds^ = dx^ + dy^ = ed^^ + 2fd^drj + gdi^^ 

e = X ^2 + y ^2 

^ = X^x^ + y^y^ 

g = + y ^2 

j = (eg - 


Truncation error of centered finite-difference approximation 


1 

12 j 


4f(Zj„,A^2 + Zj,,,Ar,2] 


Choose boundary-fitting map to 

(1) Minimize (f/j) 2 

(2) Reduce errors in separate 
I and 17 directions. 

One-dimensional monitor function* methods are 
satisfactory for (2). Boundary adjustment of 
map used for ( 1 ). 


*A. B. White, Jr., SIAM J Num Anal 16 (1979) 

C. M. Ablow and S. Schechter, J. Comp Physics 27 (1978) 
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One-dimensional monitor function takes equally 
spaced values at mesh modes. 

Simplest, effective, problem-dependent monitor 
is distance on solution surface: 

— is to be constant as i varies 

- 0 

Same for 17 direction 

< 2 > = 0 

( 1 ) and ( 2 ) plus given differential equation i*) 
for z determine solution when corner points 
0,A,B,C have been chosen. Corners are moved 

to minimize 


Examples show grids found from (1) and (2) 
for various given functions z having regions 
of sharp variation (boundary layers). In 
practice, function z and the grid mapping 
would be found by simultaneous solution of 
the complete set of differential equations 
(1), (2), and (*). 























